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Abstract. The empirical relation between the X-ray luminosity (in the 2-10 keV band) and the rate of spin-down 
energy loss L s( j of a sample of 39 pulsars is re-examined considering recent data from ASCA, RXTE, BeppoSAX, 
Chandra, and XMM-Newton and including statistical and systematic errors. The data show a significant scatter 
around an average correlation between £ x ,(2-io) an d ^sd- By fitting a dependence of i x ,(2-io) on the period P 
and period derivative P of the type L x ,(2-io) tx P a P b , we obtain a — —4.00 and b — +1.34 (i.e. a ~ —3b). This 
translates into the relation £ x ,(2-io) — £x,nor (isd/erg s -1 ) with a normalization L x ,nor = 10 -15 ' 3 erg s _1 . 
However, the reduced \ 2 is large (= 7.2) making the fit unacceptable on statistical ground. All the X-ray lumi- 
nosities lie below a critical line L x ,crit : the corresponding efficiency of conversion of rotational energy into 2-10 
keVX -rays is n x — (-^ x ,crit/-^sd) L®^^ and varies, within the sample, between 0.1 and 80%. The large dispersion 
of L x below L Xj crit indicates that other physical parameters uncorrelated with P and P need to be included to 
account for the observed emission at X-ray energies. We indicate a few possibilities that all conspire to reduce 

i X ,(2-10) • 
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1. Introduction 

About 40 of the nearly 1,400 known radio pulsars 
(Manchester 2001) have been so far detected in the X-ray 
range. Despite the great differences in their rotational pa- 
rameters, magnetic field and age, these rotation-powered 
X-ray sources seem to obey a simple empirical law. The 
data reveal a correlation between their high-energy lumi- 
nosity and their spin-down luminosity L sc j = Air 2 IP /P 3 , 
where P and P are the spin period and its derivative, and 
I is the neutron star momentum of inertia that we assume 
equal to 10 45 g cm . 

This correlation was first noticed in a small set of radio 
pulsars by Seward and Wang (1988), and later investigated 
by Becker and Trimmer (1997; BT97 hereafter) using a 
sample of 27 powered pulsars detected with ROSAT in 
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the soft energy band (0.1-2.4 keV). BT97 described the 
available data with a simple scaling relation 

L x ~ 10- 3 L sd . (1) 

However, a harder energy interval seems preferred to ex- 
plore the nature of the X-ray emission resulting from the 
rotational energy loss: at energies above ^2 keV the con- 
tribution from the neutron star cooling and the spectral 
fitting uncertainties due to interstellar absorption are re- 
duced. Saito (1998) examined the correlation between L s d 
and the pulsed X-ray luminosity using a sample of radio 
pulsars detected at high energies with ASCA (2-10 keV) 
finding the scaling relation 

£ x ,(2-io) - 10~ 21 Lf d 2 ■ (2) 

None of these studies included the statistical and sys- 
tematic uncertainties on L x resulting from the errors on 
the X-ray fluxes and on the poorly known distances of 
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most pulsars. A proper treatment of these uncertainties is 
required to quantify the strength of the correlation. 

The physical origin of the non-thermal X-ray emission 
from radio pulsars is uncertain. Polar cap models give 
predictions on the pulsed X-ray luminosity, which is at- 
tributed to inverse Compton scattering of higher order 
generation pairs on soft photons emitted by the surface 
of the neutron star and/or by hot polar caps (Zhang & 
Harding 2000). The soft tail in the inverse Compton scat- 
tering spectrum can explain the non-thermal X-ray com- 
ponent observed in many pulsars. 

Outer gap models attribute the pulsed non-thermal 
X-ray emission to synchrotron radiation of downward cas- 
cades from the outer gap particles, and include (as in the 
inner gap models) a thermal component from the hot po- 
lar caps heated by impinging particles. Both models are 
consistent with the scaling relation L x oc L s( j (Cheng, Gil, 
& Zhang 1998) in the ROSAT domain, but are in conflict 
with the observations in the ASCA (Zhang & Harding 
2000) domain where they over or under-predict the val- 
ues of L x (relative to relation [2]). In polar cap models, 
these discrepancies can be accounted for by invoking ad 
hoc model parameters for the single sources, due to the dif- 
ficulty in tracking the full cascade of particles and photons 
down to X-ray energies, coupled with the uncertainties in 
the temperature of the polar cap and of the whole surface. 

The correlation between L x and L s( j is even more diffi- 
cult to quantify when considering that pulsed emission is 
observed only in 15 of the sources. Thus, for most of the 
pulsars, the total luminosity L x can be considered only as 
an upper limit on the pulsed component. L x often con- 
tains also the (unpulsed) contribution resulting from pul- 
sar wind nebular emission (Chevalier 2000) which is in end 
related to the spin parameters P and P. 

For a comparison of the theoretical predictions with 
the observations, a statistical analysis of the data seems 
useful. Given the improved capabilities of the Chandra 
and XMM-Newton observatories, it is now timely to re- 
consider the correlation L x versus L s( j in statistical terms, 
and this is the aim of our paper. 

In §2 we describe the current sample and the crite- 
ria used to account for the uncertainties in L x . In §3 we 
present the results. In §4 we comment on the possible ori- 
gin of the large scatter seen in the data and in §5 we 
summarize our conclusions. 

2. Criteria for the data analysis 

Our aim is to estimate the fraction of the rotational en- 
ergy loss of a neutron star going into X-rays. This excludes 
from the sample the X-ray pulsars powered by accretion 
from binary companions (see e.g. Bildsten et al. 1997), 
and the Anomalous X-ray Pulsars (Mereghetti 2002), the 
luminosities of which are inconsistent with rotational en- 
ergy loss in case of neutron stars. According to these cri- 
teria, the current sample consists of the 41 pulsars listed 
in Table 1. 



Even in the case of spin-down powered neutron stars, 
some contribution(s) to the X-ray flux can derive from 
cooling. Thus we include only the following components 
in the calculation of the X-ray luminosity: (i) the non- 
thermal (pulsed) emission that originates in the neutron 
star magnetosphere; (ii) the thermal emission resulting 
from the re-heating of (parts of) the neutron star surface 
by backward accelerated particles; (Hi) the flux from the 
extended synchrotron nebulae powered by the relativistic 
particles and/or magnetic fields ejected by the neutron 
star. As a rule, when possible we subtracted the contribu- 
tion to the X-ray emission from the extended supernova 
shells. At energies ^ 2 keV, the contribution from the neu- 
tron star cooling becomes unimportant. Thus, by consider- 
ing the energy range 2-10 keV we reduce the uncertainties 
related to the subtraction of the cooling components. 

We reviewed all the data in the literature for each 
source of our sample. When good observations were avail- 
able in the 2-10 keV band, we adopted the reported flux. 
When only fluxes in a different energy range were avail- 
able, we converted them to our range, adopting the param- 
eters of the best spectral fit reported in the literature (see 
the references in Table 1). The luminosities L x ,(2-io) are 
calculated assuming isotropic emission and the distances 
listed in Table 1. 

As a main improvement with respect to previous sim- 
ilar studies, we derive a weighted fit of the L x versus L s d 
data. For this purpose, we have taken into account the 
uncertainty associated to each value of L x . Contributions 
to such uncertainties have different origins: 

(a) Uncertainties in the distances; for most of the 
pulsars, we have used the distances derived from the 
dispersion measure (DM; Pulsar Catalogue http: 
//pulsar .princeton. edu/pulsar/catalog). They 
rely on a model for the electron distribution of the 
interstellar medium (Taylor & Cordes 1993). For indi- 
vidual sources this can lead to an error in the distance 
up to a factor £ 3 (e.g. the case of PSR Jl 119-6127, 
Camilo et al. 2000, Crawford et al. 2001), but, when 
averaged over the pulsar population, the typical errors 
reduce significantly (Taylor & Cordes 1993). When 
only the distance inferred from the DM was available, 
we conservatively adopted an uncertainty of ±40%. 
This translates in an error ~ 0.3 in logL x . When 
other determinations of d exist (e.g. from parallax 
measurements or from associations with supernova 
remnants) we opted for the most accurate ones, 
assuming the uncertainties reported for the related 
measure. 

(b) Statistical errors in the number of photons; the 

errors due to photon counting statistics dominate the 
weak pulsars first detected with ROSAT, yielding un- 
certainties up to ~ 0.5 in logL x ; statistical errors in 
the count rate are negligible for brighter objects. 

(c) Errors depending on interstellar absorption; the 
conversion from the observed flux to that emitted at 
the source depends on the amount of interstellar ab- 
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Table 1. Parameters and emission of the 41 sample pulsars 



Pulsar 


P 


P 


N H 


distance 


a or kT 


Obs Band 


Detector 


y(2-10) 




ref 


PSR 


ms 


10~ 15 S8" 1 


10 21 cm -2 


kpc 


kcV 


koV 




erg s~ 1 cm~ 2 


erg s 1 




J0030+0451 


4.87 


1.00 ■ 10~ 5 


2.15 ± 0.85 


0.230 ± 0.092 


2 ± 0.2 


0.1-2.4 


PSPC 


(i.27r -) 


• io- 13 


29.88l- 73 


1 


J0218+4232 


2.32 


7.50 ■ 10" 5 


2 ± 2 


5.70 ± 2.28 


0.94 ± 0.22 


2.0-10.0 


MECS 


(4.30l° ; 2 l) 


• io- 13 


33.20+55' 32 


2 


J0437-4715 


5.76 


1.86 ■ 10" 5 




0.178 ± 0.026 


2.35 ± 0.35 


0.1-2.4 


HR.I 


(4.30ii) ' 


io- 13 


30.l9i°; 24 


3,4 


J0751+1807 


3.48 


8.00 ■ 10~ 6 


4 4 + 4 6 


2 ± 0.8 


2 ± 0.2 


0.1-2.4 


PSPC 




• io- 14 


3i-29iS: 5 8 2 


3 


J1012+5307 


5.26 


1.46 ■ 10" 5 


0.06 ± 0.01 


0.520 ± 0.208 


2.3 ± 0.2 


0.1-2.4 


PSPC 


(1.25i° S 6 ?) 


• io- 14 


29-58l°; 5 7 3 


3 


J 1024-0719 


5.16 


2.99 ■ 10~ 6 


0.2 ± 0.05 


0.350 ± 0.140 


2 ± 0.2 


0.1-2.4 


HR.I 


(8.86i 4 3 ;™) 


• io~ 15 


29.09+Q-™ 


3 


J 1744- 11 34 


4.07 


7.13 ■ 10~ 6 


0.1 ± 0.05 


n QC7+0.043 

u.oo(_ 035 


2 ± 0.2 


0.1-2.4 


HRI 


(6.44± 4 ;^) 


• io- 15 


28-97i° -; 32 


3 


B1821-24 


3.05 


1.62 ■ 10" 3 


2.9 ± 2.3 


5.1 ± 0.5 


1.89 ± 0.21 


0.7-10.0 


CIS 


(i-25lS;S) 


• IO" 12 


33.56-toH 


5 


B1937+21 


1.56 


1.06 ■ 10~ 4 


21 ± 5 


3.60 ± 1.44 


i 71 +0.05 
i- ' -0.08 


0.5-10.0 


LECS+MECS 


(3.70i° ,f ) 


• io- 13 


32.73+!!' 3 ? 

— 0.55 


6 


J2124-3358 


4.93 


1.30 ■ 10~ s 


0.35 ± 0.15 


0.25 ± 0.10 


2 ± 0.2 


0.1-2.4 


HRI 


(8-26± 3 f 8 ) 


• io- 14 


29.77l°o. 32 


7 


B0950+08 


253.07 


0.229 




0.127 ± 0.013 




0.1-2.4 


PSPC 


(2-3i U -\) ■ 


IO" 14 


28.62^0^2 


3 


B1929+10 


226.52 


1.16 


0.1 ± 0.05 


0.25 ± 0.08 


0.44 ± 0.046 c 


0.5-5.0 


SIS 


(5-6±M) ■ 


10 -14 


29.60;° ; 34 


8 


B0823+26 


530.66 


1.71 




0.380 ± 0.152 




0.1-2.4 


PSPC 


(o. 6 ;° 2 2 ) ■ 


1Q -14 


28 QQ+ 42 
^°- yy -0.62 


3 


B0114+58 


101.44 


5.85 


2.57 ± 0.2 


2.14 ± 0.856 


2.1 ± 0.2 


0.1-2.4 


PSPC 


(4.25^1») 


• IO" 15 


30.34+^;^ 


9 


B0355+54 


156.38 


4.40 


0.2 ± 0.2 


2.10 ± 0.84 


2 ± 0.5 


0.1-2.4 


PSPC 


(i-i6i 3 : 9 4 4) 


• io- 13 


31.76t?® 7 


10 


J0538+2817 


143.16 


3.67 


0.6 ± 0.6 


1.5 ± 0.6 


1.5 ± 0.5 


0.1-2.4 


PSPC 


(8.00t 8 ™) 


■ io- 16 


29.3li°; ?? 


11 


B0633+17 


237.09 


11.0 


0.13 ± 0.13 


u,1J *-0 .034 


2.19 ± 0.35 


0.7-5.0 


CIS 


(7.94+. 3 2 ) 


10 -14 


29 c.o + 0-42 
°°-0.36 


12 


B0656+14 


384.89 


55.0 


0.17 ± 0.17 


n 9s + 0.20 
°' 28 -0.10 


1.5 ± 1.1 


1.0-5.0 


SIS 


(2.05l ; 7 7 2 ) 


• io- 13 


30.26± H 


13 


B1055-52 


197.11 


5.83 


0.26 ± 0.06 


0.5 ± 0.2 


1.5 ± 0.3 


2.0-10.0 


CIS 


(1 06+ 010 ) 


• io- 14 


29-48i°; 3 | 


14 


B1951+32 


39.53 


5.84 


3.4 ± 0.5 


2.5 ± 0.2 


2.1 ± 0.3 


0.1-2.4 


PSPC 


(2-04i°- 7 |) 


• io- 12 


33.16l°; 22 


15 


B0833-45 


89.33 


1.25 • 10 2 


0.4 ± 0.1 


0.25 ± 0.03 


2.2 ± 0.4 


0.2-8.0 


ACIS-S 


(1.03± ^) 


•IO" 11 


31.86±r44 


16 


B1046-58 


123.67 


96.3 


5±1 


2.98 ± 1.19 


2 ± 0.2 


0.4-10.0 


SIS 


(2-50t° 6 5 6 8 ) 


• io- 13 


32 40+ ' 39 


17 


J1105-6107 


63.19 


15.8 


8.55 ± 5.25 


7.0 ± 2.8 


1.8 ± 0.4 


2.0-10.0 


CIS 


v — 1 . 04 ' 


• io- 13 


,Q t;r; + 0.37 

aa.oo_ 52 


18 


J 1420-6048 


68.18 


83.2 


22 ± 7 


2.0 ± 0.8 


1.6 ± 0.4 b 


2.0-10.0 


CIS 


(4-70l° ; 7 I) 


• io- 12 


oo QQ~h~0.36 
OO.OO q ^2 


19 


B 1706-44 


102.46 


93.0 


3.45 ± 3.45 


1.80 ± 0.72 


1.9 ± 0.9 


2.0-10.0 


SIS+GIS 


(1.03_ Q 24 ) 


in- 12 




20 


B1757-24 


124.90 


1.28 • 10 2 


35 ± 12 




1.6 ± 0.6 


2.0-10.0 


ACIS-S 


(7 9+°' 6 ~l • 


1Q -13 


,„ o 7 +0.20 
-0.10 


21 


J_> _L O U \-J ^ _L 


133.63 


1.34 • 10 2 


13 ± 1 


k on 4- 9 19 


2 ± 0.2 


0.1-2.4 


PSPC 


It To + 0.70\ 

(1-78I 59 ) 


i n- 13 
• 10 


qo 7 c: + 0.45 
3Z - ,O -0.70 


22 


J1811-1926 


64.67 


44.0 


13.8 ± 0.8 


7.8 ± 2.5 


1.89 ± 0.25 


4.0-10.0 


MECS 


(i.23±s-;d 


• io- 11 


34 93+ - 27 


23 


DloZo—LO 


101 .45 


75.5 


40 i 25 


A i o _l_ i an 
4.1z ± l.uo 


2 i 0.2 


0.5—2.4 


PSPC 


(1.70tt1) 


IO" 11 


*3 A r i + . 8 5 

34.51_ 1 20 


24 


B1853+01 


267.40 


2.08 • 10 2 


2.57 ± 0.2 


3.2 ± 1.3 


2.3 ± 1.1 


0.4-2.0 


CIS 


/ i r-,+ 3 \ 

(l-2Io.l) ' 


10 


33 14+ 039 

JJ ' -0.57 


25 


J2229+6114 


51.62 


78.0 


6.3 ± 1.3 


3 ± 1 


1.51 ± 0.14 


2.0-10.0 


ACIS-I 


, (in 

(i-SO-oos) 


• 10 


. . If) OQ 

33.12+g;|| 


26 


B2334+61 


495.28 


1.92 • 10 2 


2 ± 1 


2.5 ± 1 


2 ± 0.2 


0.1-2.4 


PSPC 


(4.05+. 28 ) 


10 -14 


31.46t°| 2 


27 


J0205+6449 


65.68 


1.93 • 10 2 


3 ± 2 


2.6 ± 0.6 


1.9 ± 0.2 b 


0.8-10.0 


HRC 


(l-5i U o-.3) • 


io- 11 


*ia ™ + 14 

34.08+g 21 


28 


B0531+21 


33.52 


4.21 • 10 2 


3 ± 0.5 


2 ± 0.5 


2.108 ± 0.006 


0.3-10.0 


MOS 


(9.93±°;») 


• io- 9 


36.65;° ; 2 2° 7 


29 


J0537-6910 


16.11 


51.0 


6.9 ± 3.3 


47.3 ± 0.8 


1.6 ± 0.2 a 
2.55 ± 0.15 " 


0.2-10.0 


CIS 


(B.i3li;») 


• io- 12 


36.iii°„;i 2 


30 


B0540-69 


50.53 


4.73 • 10 2 


4.6 ± 4.6 


47.3 ± 0.8 


1.83 ± 0.13 " 
2.06 ± 0.2 6 


0.2-10.0 


ACIS-I 


(3.33i«;f 9 ) 


• io- 11 


36.93lS; 2 3 3 


31 


,11119-6127 


407.75 


4.02 • 10 3 


15 ± 15 


5 ± 3 


1 4+1-0 


0.7-5.0 


CIS 


(4.74+1?) 


IO" 13 


33.13+;™ 


32 


J1124-5916 


135.31 


7.45 • 10 2 


3.17 ± 0.15 


4.8 ± 1.6 


1.9 ± 0.2 


2.0-8.0 


ACIS-S 


(i-i± : 2 2 ) • 


io- 11 


34 48+ ' 18 


33,34 


B1509-58 


150.66 


1.54 • 10 3 


12. 7± 12.7 


4.2 ± 0.5 


1.358 ± 0.014 a 
2.2 ± 0.005 " 


2.0-250 


PCA 


(l-05l° ° 3 8 3) 


• io- 10 


35.32±!!; 2 2 


35 


J1617-5055 


69.36 


1.37 • 10 2 


6.8 ± 6.8 


4.5 ± 0.9 


1.6 ± 0.3 


3.5-10.0 


CIS 


(8.86tg;«) 


• io- 12 


34.3ltg;i» 


36 


J 1846-0258 


323.60 


7.10 • 10 3 


47 ±8 


19 ± 5 


2.2 ± 0.1 


3.0-20.0 


PCA 


(3.90;° 4 ) 


IO" 11 


36.221S; 23 


37,38 



R.efs: Note : a Pulsed component photon index; b compact nebula photon index; c polar cap black body temperature. Detector : PSPC and HRI 
are instruments aboard of ROSAT; GIS and SIS are intruments aboard of ASCA; LECS, MECS and PCA are instruments aboard of Beppo-SAX; 
ACIS-I, ACIS-S and HRC are instruments aboard of Chandra; MOS is an instrument aboard of Newton-XMM. Obs Band : is the energy band of 
the observation from which the fluxes reported in column 9 arc measured or extrapolated. References : [1] Becker ct al. 2000; [2] Mineo et al. 2000; 
[3] BT99; [4] Kawai et al. 1998; [5] Saito et al. 1997; [6] Nicastro et al. 2002; [7] Sakurai et al. 2001; [8] Wang, Halpern 1997; [9] Slane 1995; [10] 
Slane 1994; [11] Sun et al. 1995; [12] Halpern, Wang 1997; [13] Grciveldinger et al. 1996; [14] Cheng, Zhang 1999; [15] Chang, Ho 1997; [16] Pavlov 
et al. 2001; [17] Pivovaroff et al. 2000; [18] Gotthelf, Kaspi 1998; [19] Roberts et al. 2001; [20] Finley et al. 1998; [21] Kaspi et al. 2001; [22] Finley, 
Ogelman 1994; [23] Torii et al. 1999; [24] Finley ct al. 1996; [25] Harrus et al. 1996; [26] Halpern ct al. 2001; [27] BT97; [28] Murray et al. 2002; [29] 
Willingale ct al. 2001; [30] Marshall, Gotthelf et al. 1998; [31] Kaaret et al. 2000; [32] Pivovaroff ct al. 2001; [33] Hughes et al. 2001; [34] Camilo et 
al. 2002; [35] Marsden et al. 1997; [36] Torii et al. 1998; [37] Gotthelf et al. 2000; [38] Mereghetti et al. 2002 



sorption, which is typically deduced (in a model depen- 
dent way) from the X-ray spectral fits. We evaluated 
the errors deriving from the poor knowledge of N# on 
a case by case basis, 
(d) Errors depending on spectral and spatial mod- 
eling; collecting area, integration time, spectral and 
spatial resolution span a large range of values for the 
observations of different sources, implying that they 
have been studied with different degree of details. This 



introduces strong differences in the uncertainties re- 
lated to the modeling of each source. So we have stud- 
ied them individually. 



In the computation of the error bars we first varied 
./Vh and the spectral parameters within the uncertainties 
reported in Table 1, taking the resulting maximum and 
minimum fluxes. The errors on fx are then propagated 
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Fig. 1. The X-ray lumi- 
nosity in the band (2-10) 
keV against the spin-down 
luminosity for the 41 sources 
of our sample. The objects 
are grouped in 5 classes 
and labeled accordingly: 
Millisecond Pulsar (filled 
squares), old pulsar (empty 
squares), Geminga-like (filled 
triangles), Vela-like (empty 
triangles) and Crab-like 
(filled stars). The solid line 
represents our equation [3]; 
it has been calculated on 
39 objects (thus excluding 
PSR J01 14+58 and PSR 
J0538+2817, also labeled 
with a circle filled with a 
cross and without their huge 
error bar; see Table 1). 



accounting for the distance uncertainties and for the error 
due to photon counting statistics. 

3. The current sample 

The 41 pulsars listed in Table 1 can be naturally grouped 
into five classes, based on their intrinsic and observational 
characteristics. 

I Millisecond Pulsars (MSPs): Six among the ten MSPs 
in our sample have been detected only by ROSAT, and 
four of these with a very low photon counting statistics 
(< 50 photons). This does not allow to asses which is 
their best spectral model. The objects studied in more 
detail with ROSAT indicate that non-thermal fits arc 
preferred (Becker & Triimper 1999). However, recent 
observations with Chandra of the MSP population in 
the globular cluster 47 Tucanae (Grindlay et al. 2001), 
indicate that a thermal component may be present at 
energies Si 1.5 keV. Considering that the possible ther- 
mal component would contribute only slightly to the 
luminosity in our energy range (2-10 keV), we have 
adopted for all the MSPs a power law spectrum with 
photon index a p h — 2 (so that the specific energy flux 
f v oc v^ a+1 ^). The list of sources belonging to this 
group is reported in the upper section of Table 1. 
II Old Pulsars: Only three of the non recycled pulsars 
in our sample have characteristics ages (t = P/2P) 
greater than 10 6 yr: PSR B0823+26, PSR B0950+08 
and PSR B1929+10. ASCA observations of the lat- 
ter two sources were reported by Wang & Halpcrn 
(1997), but it was later found that the results for 
PSR B0950+08 were affected by the presence of a 



nearby AGN (see note in Wang et al. 1998). The spec- 
trum of PSR B1929+10 in the ASCA band can be 
described by thermal emission from a small fraction 
of the star surface (^ 10 — 30 m in size). This is gen- 
erally interpreted in terms of re-heating of the polar 
caps by backward accelerated particles in the magneto- 
sphere. Assuming that the same interpretation is valid 
for PSR B0823+26 and PSR B0950+08, we use for 
these objects the flux obtained by scaling their ROSAT 
count rates to that of PSR B1929+10 and adopting 
large error bars. 

III Geminga-like: This class contains the middle aged pul- 
sars (t ~ 10 5 yr) for which the internal cooling gives 
a substantial contribution, at least at energies below 
~2 keV. For the three brightest objects of this group 
(Geminga, PSR B0656+14 and PSR B1055-52) ac- 
curate spectral modeling exist, allowing to estimate 
with a good accuracy and subtract the contribu- 
tion from the cooling to L x . For PSR B1951+32 our 
adopted L x derives mainly from the synchrotron neb- 
ula surrounding the pulsar. The remaining two sources 
(PSR B0114+58, PSR J0538+2817) are considerably 
weak preventing a detailed spectral analysis; therefore 
these two sources are not included in the statistical 
fit. However, in order to include them in the graphical 
representations of Figures 1 and 4, we have assumed a 
power-law model, with a as given in Table 1, and a ra- 
tio ~ 10 -3 between the flux in the power-law and the 
cooling component, i.e., an average between Geminga 
and PSR B1055— 52. This could largely underestimate 
the flux if it is entirely of magnctospheric origin. 

IV Vela-like: This is a rather inhomogeneous group of 
~ 10 4 — 10 5 years old pulsars. Some of them are associ- 
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Fig. 2. The X-ray luminosity in the (2-10) keV band plotted against the spin-down luminosity for 4 subclasses of objects 
in our sample. The solid lines represent the best fits of a relation of the type log£ x .(2-10) = TO l°g^sd — q calculated 
for the sources of each subclass: they are log i x ,(2-io) = (1-38 ± 0.10) log£ s d — 16.36 ± 3.64 for the Millisecond Pulsars 
(panel a), logL X: (2-io) = (1.71 ± 0.17) log L sd — 29.35 ± 6.15 for the Geminga-like sources (panel b), logL X: ( 2 -io) = 
(0.46±0.23)logL sd + 16.61±8.33 for the Vela-like sources (panel c) and logL Xj ( 2 -io) = (1.02±0.10) logLsd -3T4±3.86 
for the Crab- like sources (panel d). The dashed line is the fit given in equation [3]: see text for details. 



ated to supernova remnants. For most of these objects, 
the main contribution to £ x .(2-io) arises from their 
synchrotron nebulae, that typically have power-law 
spectra with a ph in the range 1.5-2.3 (PSR B0833-45, 
PSR B1046-58, PSR J1105-6107, PSR B1706-44, 
PSR J1811-1926, PSR B1823-13, PSR B1853+01, 
PSR J2229+61). In the case of PSR B1757-24 the 
bulk of the emission in the (2-10 keV) band seems 
of magnetospheric origin, whereas in the cases of 
PSR J1420-6048, PSR B1800-21, PSR B2334+61 the 
number of collected photons is too small for a detailed 
analysis. 

V Crab-like: These are the youngest pulsars in our sam- 
ple. They have been directly seen in the hard energy 
band, and are usually bright enough for a careful ex- 
amination of their spectra. Their emission consists of 
two components, a pulsed one of magnetospheric ori- 



gin, and a second one (unpulsed) from the synchrotron 
nebula. The pulsars of this class are listed in the lower 
section of Table 1. Only for PSR J1119-6127 and 
PSR Jl 124— 5916 observation prevents the possibility 
of disentangling the pulsed from the unpulsed com- 
ponents: for the former source due to the low photon 
statistics of the ASCA detection, for the latter object 
due to the poor time resolution of the adopted Chandra 
detector. 



4. Results 

To facilitate the comparison with previous works, we have 
first fitted 39 points in Figure 1 with a linear relation, 
obtaining 

Iog£ x ,(2-io) = 1-34 logL sd - 15.34 (3) 
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Fig. 3. X-ray luminosities in the (2-10) keV band plotted versus the rotational period (panel a), the period derivative 
(panel 6), the characteristic age (panel c), and the surface magnetic field (panel d). 



with a x 2 — 259.4 (solid line in Figure 1). The formal 
lcr intervals of uncertainty are 1.34 ± 0.03 (for the slope) 
and 14.36 ± 1.11 (for the constant term). Thus the slope 
m of the relation (||) is slightly flatter than the value of 
1.5 derived in the same energy interval for a sample of 15 
objects by Saito (1998), who considered only the pulsed 
emission. At lower energies, ROSAT data instead suggest 
a slope of TO(o.i_2.4) = 1-04 ± 0.09 for 27 sources (BT97), 
though a steeper dependence mm 1-2.4) — 1-35 was de- 
rived by Ogelman (1995) using a sample of 7 sources. At 
the intermediate energy range 0.2—4 keV Seward & Wang 
(1988) found a value of 1.39 based on a small sample of 8 
pulsars. We note however that the empirical relation of L x 
with L s d is not fully satisfied as the fit of equation (|^) is 
statistically unacceptable having an extremely large value 
of the reduced x 2 ~ 7.0. The scatter around this relation 
is remarkable, even when splitting the data in the different 
classes of sources, as shown in Figure 2. 

In the caption of Figure 2, we have reported the re- 
sults of the linear fit for 4 distinct subclasses of objects in 
our sample. We note that the Crab-like sources (panel d) 



are evenly distributed on both sides of the best-fit line of 
equation (||) (dashed line in all the panels of Figure 2), but 
their local best fit line (with a slope of m cra b = 1.02) is 
flatter than that derived from fitting all data in our sam- 
ple. Remarkably, the two sources located in the Magellanic 
Clouds (whose relative X-ray luminosity is not affected by 
uncertainties in the distance) even show an anticorrelation 
between L sd and £ x ,(2-io) : PSR J0537— 69 (whose L sd is 
about 3 times greater than that of PSR B0540— 69) ap- 
pears about ten times dimmer in the ASCA band than 
PSR B0540— 69: we note that the ratio between the two 
luminosities, for these two objects, scales with the ratio of 
their period derivatives P. 

Geminga and the other cooling pulsars have a rela- 
tively steep slope m gem = 1.71 and (with the exception of 
PSR B0355+54) are underluminous of 1 — 1.5 dex relative 
to the dashed line (eq. [3]) describing the whole sample. 
This could suggest that either the (subtracted) thermal 
component has been overestimated or that these sources 
are preferred targets for the mechanisms which could re- 
duce the luminosity in the 2-10 keV band (see §5). 
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Fig. 4. The X-ray luminos- 
ity in the band (2-10) keV 
against the spin-down lumi- 
nosity for the 41 sources of 
our sample. Objects and la- 
bels are as in Figure 1. The 
dot-dashed line represents the 
line logL Xi ( 2 _io) = log£ s d, 
whereas the solid line is the 
critical line of equation [6], 
obtained excluding the two 
pulsars PSR J0114+58 and 
PSR J0538+2817. All points 
in this diagram locate below 
the solid line. 



The typical spin down ages, the morphologies and the 
spectral characteristics of the Vela-like pulsars classify 
them between the Crab-like and the Geminga-like objects. 
Thus, it is not surprising that the data plotted for this 
group of sources (panel c of Figure 2) display features of 
both the Crab-like and the Geminga-like plots. In fact the 
slope TOyela = 0.46 is much flatter than that of the entire 
sample and a significant number of these objects appear 
much dimmer than predicted by equation [| 

Finally, Figure 2 shows that the MSPs are the sources 
which better follow the general logi x ( 2 -io) vs logi s( j re- 
lation, both in terms of slope m mS p = 1-38 and of lumi- 
nosities. Interestingly, the MSPs in 47 Tuc observed in 
the (0.5-2.5) keV interval with Chandra appear to have a 
weaker dependence of L x on L s d than those in the field, 
with a best median estimate for the slope of 0.55 ± 0.2 
(Grindlay et al. 2002). 

Figures 3a and 3b illustrate that the scatter in the 
values of £ x ,(2-io) (already appearing in Figure 1) is even 
larger when the X-ray luminosities are plotted separately 
versus the periods and period derivatives. We also note 
the lack of any correlation between the plotted quantities. 
There are two physical ways of combining P and P to 
give the characteristic age (r = P/2P) and the surface 

magnetic field (B s = 3.2 x 10 19 VPP G) of a pulsar. In 
Figures 3c and 3d we report £ x ,(2-io) versus these two 
quantities: while the spread in the data remains large, a 
correlation (already noted by BT97) is visible in £ x ,(2-io) 
versus r, when the subsample of the MSPs is excluded. 

This suggests the possibility to improve the quality of 
the fit by using a more general function of P and P. In 
fact, the proposed physical models for the X-ray emission 
rely on mechanisms depending on different combinations 



of these two quantities. Therefore, we have explored a fit 
of the type 

log£ x ,(2-io) = alogP + MogP + c . (4) 

Using a x 2 minimization code, we found a = — 4.00±0.13, 
b = 1.34 ± 0.03, c = 47.11 ± 0.32 with a X 2 = 259.3, 
implying that the fit is still unacceptable. Only enlarg- 
ing the error bars on logL Xj (2_io) up to absolutely un- 
reliable values of ~ 0.9 (in units of log(erg s -1 )), this 
modeling of the data would become statistically accept- 
able (with x 2 /d.o.f. = 1.1). Rearranging the best-fit, 
we get logL X)(2 _io) = 1.34(-2.991ogP + logP + c') = 
1.34(-3 log P + log P + c') + 0.02 log P. Since logP spans 
the interval — 3 — > 1, the last addendum is of the order of 
the la uncertainty on the constant term, and thus we can 
approximate the formula, getting 

log£ X) (2-io) - 1.34 log L sd - 15.30 (5) 

which is very close to the best fit of equation (|J). This 
exercise has shown that, when properly accounting for the 
uncertainties in the measured £ x ,(2-io)j 110 combination of 
P and P of the type P a P b can fit the data in a statistically 
acceptable way. However, it is remarkable that, among 
all the possible combinations of P and P, the one which 
better describes the data selects again a scaling with L s d- 
We note that all the data of Figure 4 lie under the 
critical line 

L XjCrit = 10- 18 - 5 (-^4) ergs" 1 (6) 
\erg s _i / 

obtained searching for that line having the minimum 
weighted distance to the observed points. This line gives 



8 



Possenti et al.: Re-examining the L x versus L s d correlation of Pulsars 



0) 



O 



M 
O 



CO 
CO 



CO 
CO 



CO 



C\2 
CO 



o 

CO 




34 



□ 



□ 



□ 



□ 



□ 



□ 



□ (2-10 keV) Luminosity 
* > 400 Mev Luminosity 



35 



36 37 
Log L gd (erg s _1 ) 



38 



39 



Fig. 5. High energy lumi- 
nosity as a function of L s d 
in two selected energy do- 
mains: in the (2-10) keV 
(L( 2 -io)kev) an d above 400 
MeV (X ( >4oomcV))- Stars rep- 
resent young 7-ray pulsars 
taken from Zhang & Harding 
(2000; Fig. 3); open squares 
refer to the pulsars of our 
sample (with the exclusion of 
MSPs). The dot-dashed line 
gives the locus where lumi- 
nosity equals the spin-down 
power. The dashed lines gives 
the critical luminosities (de- 
scribed as single power laws) 
in the two energy domains 
while the broken solid line in- 
terpolates the two lines, and 
indicates the observed change 
in slope, as discussed in sec- 
tion 5. 



the maximum efficiency of conversion of spin-down energy 
in X-rays: 

T / T \ 48 

^t =10 -18.5/_^^ (?) 

L s d V erg s V 

For L s d in the interval 10 32 — 10 38 erg s _1 , the efficiency 
varies within 10~ 3 < rj x < 0.8. 

5. Discussion 

L x ,crit can be interpreted as the line giving the maximum 
luminosity that can be attributed exclusively to magneto- 
spheric processes, and the large scatter below that line as 
due to additional mechanisms that tend to reduce L x and 
that are independent of P and P. Note that the width of 
the spread in the values of L x does not vary over the whole 
interval of L s d, thus supporting further this interpretation. 

The magnetospheric emission is pulsed and intrinsi- 
cally anisotropic. Thus, the detection of a pulsed sig- 
nal is susceptible to various geometrical corrections. In 
the computation of L x we have assumed that the entire 
beam is seen. The angle a between the spin axis and the 
dipolar magnetic field controls the intrinsic width of the 
emission cone. According to Zhang & Harding (2000), 
Z/ x does depend on a. This effect is particularly pro- 
nounced for a=90° and causes a drop in the luminosity. 
Disappointingly, it is still difficult to have a direct mea- 
sure of a, which is known, with large uncertainty, only for 
a handful of sources (Miller & Hamilton 1993). The Vela 
pulsar is recognized to be an orthogonal rotator: this can 
explain the drop of L x by <~ 1 — 2 orders of magnitude 
relative to L xcrltl making the observation consistent with 
the prediction, within the error determination. However, 



the small a ~ 30° of PSR 0656+14 is not sufficient to 
explain its displacement from I/ X>cr it- It is thus plausible 
that a influences the values of L x but is not sufficient to 
explain the scatter seen in the L x versus L s d correlation. 

Environmental effects can also be a source of scatter 
for the young pulsars of the sample. As an example, PSR 
0537—69 is one order of magnitude less luminous than the 
Crab pulsar despite the similarity in the value of L s( j and 
in the apparent morphology. Wang et al (2001) explain 
the smaller efficiency in terms of a stronger confinement 
of the pulsar nebula by a surrounding star formation re- 
gion, which combined with the high pulsar velocity (600 
km s _1 ) gives a smaller shock radius of a fraction of a 
parsec. The relativistic particles responsible for the X-ray 
emission thus escape the nebula in a time much shorter 
than the time of radiative energy loss. 

The 2-10 keV energy band in which the sources are 
detected is selected by instrumental needs, but it is not 
necessarily optimal to sample the bulk of the broad-band 
high energy X-ray emission from pulsars. This is what 
emerges from a spectral analysis of a number of bright 
pulsars carried out with BeppoSAX data (Massaro et al. 
2000). These new observations indicate that the photon 
spectral index (a) of the non-thermal emission is itself a 
function of energy a(E) = a + 2b\og(E/E ) (with a and 
b parameters of the spectral fit and Eg = 1 keV) that is 
consistent with the observational data when extrapolated 
in the UV-optical-IR range and at gamma-ray energies 
below 30 MeV. The spectral photon distribution peaks at 
a characteristic energy E max = E ■ 10~ a / 4fc th from source 
to source. For Crab, E max ~ 21 keV so that the 2-10 
keV interval contains a large fraction of L x . For Vela, the 
maximum would fall at about one MeV and then it would 
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imply a low L x , just as observed (relative to L x , c rit)- The 
existence of a characteristic energy or, more generally, of 
an energy threshold could cause the deviation of L x from 

The 7-ray luminosity (above 400 MeV) seems to fol- 
low a relation of the type L 7 oc L|? d 5 (Thompson et 
al. 1997; Thompson 2001). In Figure 5 we have drawn 
(in analogy to £ x ,crit m the 2-10 keV band) the crit- 
ical 7-ray luminosity (assuming emission in one stera- 
dian) L 7 . cr it ~ 10 16 ' 6 (£ s d/erg s _1 ) 5 erg s _1 derived from 
Figure 5' of Zhang & Harding (2000). We note that at 
least three young pulsars (MSPs are excluded in this dis- 
cussion), with high L s( j lie in between L XjCr it and L 7iCr it- 
This indicates that for these sources the X-ray luminosity 
(in the interval 2-10 keV) dominates over L 7 , and that the 
correlation with i s( j of the luminosity integrated above 1 
eV is steeper than indicated by Thompson (2001) above 
a value of L s d which is around 10 36 — 10 36 5 erg s _1 (as 
indicated in Figure 5). A possible interpretation of this 
steepening is that the physical mechanisms at play favor 
the emission in the X-ray channel when L s d exceeds some 
value: cascade processes can develop fully so that higher 
generation pairs can produce photons at energies of only 
a few keV. In this framework one expects that the slope 
n of the efficiency which is a function of L s( j and of the 
selected energy band r\i, oc , (where v is the centroid 
frequency of the observational band) gradually increases 
when moving toward less energetic frequencies v. Since the 
observed efficiency in the 7-ray energy band (above 400 
MeV) is rj-y oc L~®' b while in the 2-10 keV range varies as 
?7x oc ig d 5 , we predict that the maximum efficiency should 
become almost independent of L s( j (i.e., n ~ 0) in an en- 
ergy band intermediate the two. Future observations of 
pulsars in the energy range of INTEGRAL between (0.1- 
10) MeV may help in testing this prediction. 

6. Conclusions 

The analysis of the current data (with their uncertainties) 
on the non-thermal X-ray emission from rotation-powered 
pulsars shows: 

1. No monomial combination of P and P fits £ Xl (2-io) m 
a statistically acceptable way. 

2. Still, a correlation between the X-ray luminosity in the 
band 2—10 keV and L s( j persists in the data and the 
preferred scaling is £ x ,(2-io) P — 4 - 00 P 1 



-4.00 pi. 34 _ ri.34 

^sd ■ 



3. All the data lie below a critical line L x , C rit oc L^, 
providing the maximum efficiency of conversion of L s d 
in X-ray emission. 

4. Geometrical effects, cut-off energy scales and envi- 
ronment could be the causes of the reduction in the 
detected i x ,(2-io)i thus explaining the large scatter 
present in the data. 

After comparing our results with those derived sim- 
ilarly in the 7-ray band, we suggest that the maximum 
efficiency of conversion of spin-down power in high energy 



emission should be almost independent of L s( \ in the en- 
ergy band that will be explored by INTEGRAL. 
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